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PREFACE 


This  Memorandum,  prepared  for  the  Advanced  Research  Projects 
Agency,  Is  part  of  a  study  of  those  phenomena  which  affect  the  per¬ 
formance  of  optical  or  infrared  reconnaissance  and  guidance  equipment. 
The  objective  of  these  studies  is  to  provide  sufficient  understanding 
for  the  system  analyst  to  compute  performance  estimates  under  various 
operational  conditions. 

A  quantitative  understanding  of  the  effect  of  atmospheric  turbu¬ 
lence  on  a  beam  of  light  of  finite  cross  section  is  required  for  the 
prediction  of  the  performance  of  various  devices  employing  lasers  for 
target  acquisition  or  guidance  in  tactical  missions.  Such  applications 
are  characterized  by  near-horizontal  propagation  paths  near  the  ground 
of  the  order  of  one  to  tens  of  kilometers  in  length.  This  Memorandum 
discusses  the  evidence  in  support  of  a  log-normal  distribution  for 
the  fluctuating  irradiance  of  light  after  propagation  through  a  tur¬ 
bulent  atmosphere,  calculates  the  modulation  transfer  function  and 
phase  structure  function  for  the  wave,  demonstrates  reasonable  agree¬ 
ment  between  the  computed  results  and  experiments,  and  points  out  the 
significance  of  these  results  in  determining  the  signal-to-noise  ratio 
of  an  optical  heterodyne  communications  system  and  in  seeing  through 
the  atmosphere. 

These  results  should  be  of  use  to  those  interested  in  tactical 
applications  of  laser  range  finders,  laser  line  scanners,  and  the 
various  guidance  systems  employing  an  illuminating  beam. 


SUMylARY 


Tatarski's  analysis  of  the  first-order  (single  scattering) 
distribution  function  based  on  the  Rytov  approximation  is  shown  to 
be  in  error,  and  the  recent  analysis  of  deWolf,  which  yields  a 
Rayleigh  distribution  for  the  fluctuating  intensity  of  light  in  a 
turbulent  atmosphere  in  the  limit  of  long  optical  propagation  paths, 
is  disputed.  The  experimental  evidence  for  a  log-normal  distribution 
for  the  intensity  fluctuations  is  pointed  out. 

We  conclude  that  none  of  the  theories  proposed  can  account,  in 
a  consistent  fashion,  for  the  experimentally  observed  far-field 
irradiance  statistics. 

Expressions  are  derived  for  the  modulation  transfer  function  and 
the  phase  structure  function.  The  dependence  of  these  wave  properties 
on  the  turbulence  parameters  within  the  Kolmogorov  inertial  subrange 
is  found  and,  in  particular,  it  is  demonstrated  that  certain  optical 
measurements  of  the  phase  structure  function  and  the  modulation  trans¬ 
fer  function  can  be  interpreted  in  terms  of  an  outer  scale  of  tur¬ 
bulence  of  the  order  of  10  to  50  cm,  values  which  typically  can  be 
expected  under  inversion  conditions.  For  sufficiently  short  paths, 
it  is  pointed  out  that  the  theoretically  achievable  resolution 
through  the  atmosphere  and  the  signal-to-noise  ratio  of  an  optical 
heterodyne  detection  system  can  be  larger  than  previous  calculations 


have  Implied. 
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I.  INTRODUCTION 


The  properties  of  an  optical  wave  which  has  propagated  along  a 
path  through  the  atmosphere  are  of  considerable  interest  in  communi¬ 
cations,  laser  radar  systems,  and  general  applications  where  a  know¬ 
ledge  of  the  degradation  of  resolution  affects  the  performance  of  a 
system. 

The  probability  distribution  for  the  intensity  fluctuations  at 
a  fixed  point  in  the  turbulent  medium  is  discussed  in  the  next  section. 
No  new  theoretical  results  are  established;  however,  the  analysis  of 
Tatarski  is  shown  to  be  in  error,  the  results  of  deWolf,  yielding  a 
Rayleigh  intensity  distribution,  are  also  disputed,  and  the  experi¬ 
mental  evidence  in  support  of  a  log-normal  distribution  is  pointed 
out. 

In  the  third  section  we  consider  wave  characteristics  related 
to  the  correlation  of  amplitude  and  phase  over  distances  transverse 
to  the  direction  of  propagation.  Expressions  are  derived  for  the 
wave  phase  structure  function,  defined  as  the  mean  of  the  squared 
phase  difference  as  a  function  of  transverse  direction,  and  the 
modulation  transfer  function  (MTF) ,  which  is  the  autocorrelation 
function  of  the  complex  field  in  the  transverse  direction.  This 
latter  quantity  is  important  in  determining  the  limiting  resolution 
obtainable  in  forming  an  image  through  a  turbulent  medium,  and 
enters  into  calculations  of  the  signal-to-noise  ratio  of  an  optical 
heterodyne  system.  The  dependence  of  these  wave  properties  on  the 
turbulence  parameters  within  the  Kolmogorov  inertial  subrange  is 
obtained  and,  in  particular,  it  is  demonstrated  that  certain  optical 
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measurements  of  the  phase  structure  function  and  the  NTF  can  be 
Interpreted  in  terms  of  an  outer  scale  of  turbulence  of  the  order  of 
10  to  50  cm,  which  can  be  expected  under  inversion  conditions.  This 
result  of  the  analysis  suggests  that,  for  sufficiently  short  paths, 
the  theoretically  achievable  resolution  and  signal-to-noise  ratios 
can  be  larger  than  previous  calculations  implied. 
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II.  IRRAPIMCE  STATISTICS 

The  probability  distribution  for  the  intensity  fluctuations  of 
an  optical  wave  in  a  turbulent  medium  has  received  considerable 
attention  in  the  recent  literature. The  two  theoretical  treat¬ 
ments  of  the  problem  of  determining  the  distribution  function  which 
have  received  the  most  attention  are  those  due  to  Tatarski  and 
deWolf.  Tatarski  has  claimed  the  distribution  function  in  I  for 
plane  waves  in  the  far  field  of  the  scattering  eddies  (where  geometric 
optics  is  no  longer  applicable),  but  still  in  the  single  scattering 
zone,  to  be  log-normal  (i.e.,  £nl  is  normally  distributed).  deWolf, 
on  the  other  hand,  has  claimed  the  distribution  function  to  be  a  Rice 
distribution,  in  general,  which  approaches  a  Rayleigh  distribution  in 
the  limit  of  many  scatterings.  The  experimental  evidence,  under 
conditions  where  the  wind  speed  is  sufficiently  high  that  the  tur¬ 
bulence  is  fully  developed,  does  seem  to  bear  out  the  log-normal 
relationship  in  both  the  single  and  multiple  scattering  regimes. 

In  this  section  we  would  like  to  take  issue  with  both  of  the 
theoretical  treatments  mentioned  above,  notwithstanding  the  fact 
that  Tatarski 's  calculations  appear  to  predict  the  same  form  of  the 
distribution  as  indicated  by  experiments. 

Tatarski's  solution  is  based  on  assuming  a  solution  to  the 
scalar  wave  equation 

+  k^(l  +  nj^)^”!  U  =  0  (1) 

of  the  form  U  =  e'^,  where  k  is  the  wave  number  of  the  unperturbed 
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wave,  nj^  (ln^|  «  1)  is  the  fluctuating  part  of  the  index  of  refrac¬ 
tion,  and  U  represents  a  component  of  the  electric  or  magnetic  field. 

CO  t  h 

The  function  ^  is  then  expanded  as  where  the  i  term  i|r^  is 

proportional  to  the  i^^  power  of  n^^.  The  function  the  so-called 

first-order  Rytov  solution,  can  be  recognized  as  Uj^/Ug,  where  is 

th  ” 

the  i  Born  solution  to  Eq.  (1)  obtained  by  expanding  U  =  S  U. . 

ikz  " 

For  plane  waves ,  Ug  =  e  ,  and 


ikz' 


(2) 


Applying  the  central  limit  theorem  to  Eq.  (2),  Tatarski  concluded  that 
and  hence  ijf^,  was  normal.  The  intensity 


I  =  UU*  =  lUgl^  j^l  +  2Re(UgUp 


(3a) 


,2Reil)i- 


(3b) 


will  be  normal  if  Eq.  (3a)  is  used,  and  log-normal  if  Eq.  (3b)  is 
considered.  Obviously,  at  least  one  of  these  expressions  yields  an 
incorrect  distribution  to  order  n^^.  Yura^'^^  has  recently  shown  that 
it  is  necessary  to  retain  the  second-order  Rytov  term  in  order  to 
compute  the  standard  deviation  of  the  intensity  fluctuations  correct 
to  first  order  in  n^^.  It  follows  that  the  irradiance  statistics  and 
the  distribution  function  can  only  be  determined  to  first  order  in 
n^  in  the  Rytov  formalism  if  the  dependence  of  I  on  is  obtained. 

In  order  to  make  the  argument  more  general,  let  I  =  f(x),  where 
f(x)  depends  on  the  form  of  the  substitution  used  to  solve  the  wave 
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equation  (e.g.,  in  the  Rytov  formalism,  f(x)  =  e  )  and  expand 

CO  i 

X  =  S.  where  x.  ~  n, .  Then,  correct  to  second  order  in  n,  , 

1=0  111  1 


-Ft  \  j_  Ml  t  j.  ^ 

ox 

0 


f^x)  =  f2(xQ)  +  2f(xQ)  m  (x^  +  xp 

^0 


+  ^  ^1 
^0  Xq 


Conservation  of  energy  requires  f(x)  =  Iq  =  fC^ig)  (to  all  orders  in 
n^)  ,  and  utilizing  =  0  yields 


^1  r  = 
0 


_  1  aii  ^2 

2  Bx2  1 
^0 


Substituting  from  Eq.  (5)  yields  for  the  standard  deviation 


aj  =  |cf(x)  -  fWl^l  )(m 


One  can  then  inquire  if  there  is  any  consistent  first-order  function 


f(xQ  +  x^)  ,  where  x^^  =  0,  which  will  satisfy  energy  conservation  and 

yield  the  correct  standard  deviation  to  order  n..  ,  Eq.  (6).  Obviously, 

2  ^ 

energy  can  be  conserved  only  if  ^—r  =  0. 

dx'^ 

^0 
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Substituting  in  Eqs.  (4a)  and  (4b)  then  yields  Eq.  (6)  for  the 
dispersion.  Hence,  the  only  physically  consistent  first-order  function 
which  yields  the  correct  standard  deviation  to  order  n^  is  a  linear 
one,  corresponding  to  the  Born  approximation.  Therefore,  Tatarski's 
derivation  of  a  log-normal  distribution  for  propagation  paths  domi¬ 
nated  by  single  scattering  is  invalid.  The  correct  first-order 
statistics  are  given  only  by  the  first-order  Born  approximation, 
which  implies  a  normally  distributed  intensity  distribution.  This 
result  should  be  valid  for  propagation  paths  small  compared  to  the 
single-scattering  length,  but,  as  will  be  pointed  out  in 

further  detail,  the  experimental  evidence  supports  a  log-normal 
distribution. 

A  calculation  supposedly  valid  for  all  propagation  distances 
has  been  made  by  deWolf.^^^  Starting  with  the  high-frequency  approxi¬ 
mation  to  the  Fourier-transf ormed  wave  equation,  and  assuming  the 
refractive  index  fluctuations  to  be  a  gaussian  random  process,  deWolf 
purports  to  have  calculated  the  distribution  function  by  computing 
the  intensity  moments  of  all  orders.  With  the  gaussian  assumption, 
the  2ra-point  correlation  breaks  down  into  a  sum  of  products  of  two- 
point  correlations.  J.  P.  Laussade  claims  that  deWolf  arbitrarily 
discards  all  but  one  of  the  products  in  the  sum.  In  Laussade 's 
formulation,  the  integrals  for  the  higher  moments  could  not  be  repre¬ 
sented  in  terms  of  the  lower  order  moments,  and  no  formula  corre¬ 
sponding  to  deWolf's  Eq.  (22)  could  be  found. 

Laussade,  J.P. ,  "Theoretical  Study  of  Optical  Wave  Propagation 
Through  a  Turbulent  Medium  and  its  Applications  to  Optical  Communi¬ 
cation,”  Doctoral  Thesis,  California  Institute  of  Technology,  Pasadena, 
California,  November  1968. 
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For  strong  turbulence,  or  long  paths,  deWolf's  moments  tend 

toward  those  of  a  Rayleigh  distribution.  deWolf  has  considered  the 

function  F^^[F(^nI)]  where  F^  is  the  cumulative  distribution  for  a 

normalized  gaussian  random  variable  and  F  is  the  same  function  for 

an  arbitrary  random  variable.  Plotted  versus  log  I,  this  function 

will  be  a  straight  line  if,  and  only  if,  I  is  log-normally  dis- 

tributedo  deWolf  then  "manufactures"  data  which  are  supposedly 

Rayleigh  distributed  in  I,  and  claims  the  above  function  cannot  be 

distinguished  from  a  straight  line  over  the  four-octave  range  of 

(2) 

data  of  Gracheva.  We  do  not  know  how  he  manufactured  his  data, 

but  he  could  have  plotted 


cerf 


exp 


2v. 


2ct: 


} 


versus  v  (=  Xnl) 


which  is  the  exact  expression  for  F^^[F(v)],  and  observed  the  range 

over  which  the  curve  appeared  to  be  linear » 

Rather  than  follow  this  procedure,  we  plotted  the  cumulative 

distributions  for  I  distributed  according  to  a  log-normal,  Rayleigh, 

and  exponential  distribution  (the  latter  corresponding  to  a  Rayleigh 

amplitude  distribution)  in  Fig.  1.  The  curves  were  plotted  on  the 

(4) 

same  axis  over  a  two-decade  range  for  which  Fried,  et  al. ,  have 
taken  data.  There  is  only  one  parameter  in  the  Rayleigh  and  ex¬ 
ponential  distributions,  and  in  each  case  it  was  chosen  such  that 
I,  the  mean  intensity,  corresponded  to  the  mean  in  Fried's  data. 

Energy  conservation  requires  I  to  equal  the  intensity  in  the 


absence  of  turbulence.  The  difference  between  these  distributions 


Fig.  1--Cumulative  log-normal,  exponential,  and  Rayleigh  distributions  with  equal  means 
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for  intensities  greater  than  two  standard  deviations  from  the  mean 
is  obvious.  A  quantitative  measure  of  these  differences  could  be 
obtained  from  a  chi-squared  test  for  each  of  these  distributions, 
but  this  hardly  seems  necessary,  as  in  this  region  we  are  essentially 
comparing 


with  e 


2 


a  readily  observable  difference. 

In  conclusion,  we  feel  that  at  present  the  only  consistent 
theoretical  treatment  of  irradiance  statistics  in  the  far  field  is 
obtained  from  the  Born  approximation,  which  yields  a  normal  distri¬ 
bution.  It  should  be  pointed  out  that  geometric  optics,  valid  in 
the  near  field,  does  yield  a  log-normal  intensity  distribution. 
However,  in  addition  to  the  experiments  reported  here,  experiments 
at  Boulder  over  propagation  paths  of  up  to  145  km  bear  out  the 
log-normal  distribution.  Hence,  a  modified  geometric  optics  approach 
might  lead  to  a  better  understanding  of  this  phenomenon. 


Private  communication  from  Dr.  R.  S.  Lawrence,  Environmental 
Science  Services  Administration,  Boulder,  Colo. 
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III.  WAVE  STRUCTURE  AMD  MODULATION  TRANSFER  FUNCTIONS 


With  regard  to  our  interpretation  of  the  measurements  of  D  and 

tP 

the  MTF,  we  note  the  following  formulas;  First,  the  complex  field  at 


the  point  can  be  written,  correct  through  terms  of  second  order,  in 


the  fluctuations  of  refractive  index,  n 


(5) 


U(r.)  =  U  exp 
—1  o 


+  cp2(%)  -  I 


(7) 


where  is  the  field  in  the  absence  of  turbulence,  is  the  (first- 

^1 

order)  Rytov  approximation  (=  —  where  U.  is  the  first-order  Born 

approximation)  ,  and  =  — ,  where  is  the  second-order  Born  approxi- 

o 

mation.  Noting  that  conservation  of  the  average  energy  through  second 


order  in  n^^  implies  Re  cp2  +  =  0,  we  obtain  for  the  MTF,  to 


second  order  in  n^. 


M 


=  U(r^)  UCr^)*  =  exp  ^  li|;]^(r^)  -  1^(12)  1^]- 


(8) 


where,  for  convenience,  we  have  taken  lU^^I  =  1;  ^nd  the  asterisk 
denotes  the  complex  conjugate  of  the  appropriate  quantity.  Writing 
+  icp(t_^)  j  where  i  is  the  log-amplitude  and  cp  is  the 
phase  of  the  wave,  we  obtain 


M 


=  exp  {-  I  (D^  +  D^)} 


(9) 


where  and  are  the  log-amplitude  and  phase  structure  functions. 


A  bar  over  a  quantity  indicates  the  ensemble  average  of  that 
quantity. 
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respectively.  In  the  case  of  isotropic  homogeneous  turbulence,  Eq. 
(9)  can  be  written  as 


M  = 


exp 


{-  [C(0) 


(10) 


where  C(p)  is  the  sum  of  the  log-amplitude  and  phase  covariance 
functions,  and  p  =  ”  ^2 1 ’ 

(3) 

Tatarski  calculates  the  phase  and  log-amplitude  structure 
functions  in  a  plane  perpendicular  to  the  direction  of  propagation 
of  a  plane  wave,  with  wave  number  k  which  has  propagated  a  distance 
L,  to  be 


D  (p)  =  4n^k^L 

cp 

I 


$  (K)  KdK 
n 


(11) 


where  J  is  the  Bessel  function  of  zero  order,  and  5  (K)  is  the  three- 
o  n 

dimensional  spectral  density  of  refractive  index.  Substituting  Eq. 
(11)  into  Eq.  (9)  yields  the  MTF  for  plane  waves: 


M(p)  =  exp 


r  2  2  r 

-  -  4Tr  k  L 
^  Jr 


[1  - 


Jo<Kp)] 


(K)  KdKV 


(12) 


Equation  (12)  has  been  shown  to  be  correct  to  all  orders  in  n^^  for  a 
gaussian  random  process,  and  has  been  demonstrated  valid  through 
second  order  in  n^^  for  arbitrary  $^(K) .  For  spherical  waves,  J^(Kp) 
in  Eq.  (12)  is  replaced  by  J^(Kpx)dx. 

The  Kolmogorov  theory  yields  for  the  refractive  index  structure 
function: 


Laussade,  op.  dt. 


12 


D_^(r) 


2  2/3 

C  r 
n 


for  I  «  r  «  L 
o  o 


i'T)  ^  ^ 

n  o  \Jl  /  o 

o 


(13) 


where  £  and  L  are  the  inner  and  outer  scales  of  turbulence,  respec- 
0  0 

tively,  and  the  multiplicative  coefficients  were  chosen  to  make 

D  (r)  continuous  at  r  =  £  .  Tatarski  demonstrates  the  insensitivity 
no 

of  D,  and  D  to  the  behavior  of  $  (K)  for  K  <  1/l  when  /XL  «  L  • 

£  cp  n  o  o 

Hence  the  refractive  index  spectral  density  (assuming  isotropic 
turbulence) 


).033c2  ,-(M„/5.92)' 


#„(K)  - 


+  i, 


■') 


11/6 


(14) 


extrapolated  from  Tatarski,  which  implies  a  flat  spectrum  for 
K  <  l/i-Q  (^Q  ~  L^/2tt),  should  give  an  adequate  quantitative  depen¬ 
dence  of  D  and  M  on  the  turbulence  parameters. 

In  Figs.  2  and  3  we  have  reproduced  some  experimental  curves 

of  D  for  propagation  paths  of  0.5  and  3.5  km,  respectively.  On  the 
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same  axis  we  have  plotted  D  ,  obtained  by  numerically  integrating  Eq. 

9 

(11)  for  plane  waves  (which  correspond  to  the  experimental  situation 

in  Ref.  6).  The  results  are  very  insensitive  to  £^ ,  and  a  nominal 

value  of  0.1  cm  was  assumed.  The  phase  structure  function  will  scale 
2 

as  C^,  and  will  be  roughly  determined  by  the  transverse  distance 

corresponding  to  the  inflection  point  in  D  .  Comparison  shows  that 

2  -16  -2/3 

reasonable  agreement  can  be  obtained  with  C  -^(5X10  )cm 

n  ' 


and 


13 


Fig,  2--Coniparison  of  computed  phase  structure  function  with  data  of 
Bertolotti,  et  al , ,  at  a  range  of  0.5  km 
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~  12-30  era.  Larger  values  of  would  not  yield  a  leveling  off 

of  D  within  ~  10  cm  to  the  extent  obtained  in  the  experiment.  If 
CD 

the  measurements  had  been  made  in  the  daytime  rather  than  after 
sunset,  larger  L^'s  and  a  correspondingly  greater  range  of  validity 
for  the  "5/3  law”  would  be  expected. 

In  Fig.  4,  Djurle  and  Back's^^^  experimental  curves  of  the  MTF 
for  a  path  length  of  11  km  have  been  reproduced.  With  =  0.1  cm, 
the  curves  can  be  fairly  well  represented  by  using  spherical  waves 
and  choosing  ~  10  cm  and  25-50  cm.  The  interpretation 

of  the  shape  of  the  MTF's  implying  outer  scales  of  the  magnitude 
suggested  depends  on  the  experimental  values  being  correct  to  within 
~  10  percent.  It  is  noteworthy  that  no  estimates  of  the  experimental 
errors  are  given  for  either  D  or  the  MTF.  This  seems  to  be  character- 
istic  of  most  optical  experiments,  which  makes  theoretical  inter¬ 
pretation  difficult. 

Measurements  made  under  temperature  lapse  conditions^  ^  are 
usually  inferred  to  yield  an  outer  scale  of  turbulence  of  the  order 
of  a  meter  (when  the  path  is  ~  1  meter  above  the  ground).  However, 
most  optical  experiments,  including  the  ones  discussed  here,  are 
performed  after  dark  when  the  probability  of  temperature  inversion 
is  high.  During  the  semi-stable  conditions  occurring  during  in¬ 
version,  a  reduction  in  turbulent  energy  generation  takes  place. 

The  usual  picture  presented  is  the  cascading  of  energy  from  larger 
to  smaller  turbulent  eddies,  and  with  a  decrease  in  energy  generation 
a  reduction  of  the  outer  scale  size  would  occur.  This  phenomena  is 
supported  by  the  experiments  of  Deitz^^^  and  Tsvang.^^^  Nighttime 


Modulation  transfer  functian 


Fig.  4--Conipari son  of  computed  modulation  transfer  function  with  data 
of  Djurle  and  Back  at  a  range  of  11  km 
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measurements  of  the  temperature  structure  function  between  two  sensors 
lead  Deitz  to  infer  an  outer  scale  of  30  cm  or  less.  Tsvang's  mea¬ 
surements  display  a  shifting  of  the  center  of  gravity  of  to  larger 
K  as  the  vertical  temperature  gradient  changed  from  lapse  to  inversion 
conditions. 

One  can  argue  on  physical  grounds  that  a  lower  limit  to  the  MTF 
must  be  obtained  when  the  transverse  distance  p  is  sufficiently 
large  that  the  optical  wave  arriving  at  the  points  and  been 

scattered  through  (statistically)  independent  media.  In  this  case, 

^  Vf  _0t  /t  Vf 

as  p  00,  M  =  U(r^)  UCr^)  U(r^)  UCr^)  —  e  Hence,  if  the 

optical  path  is  sufficiently  short,  the  mediimi  does  not  limit  the 
resolution  which  can  be  obtained. 

A  decrease  in  the  inertial  subrange  could  thus  explain  the 
"saturation"  in  the  phase  structure  function  and  the  MTF.  If  the 
limiting  value  of  the  MTF  were  sufficiently  high,  optical  resolution 
could  be  improved  by  employing  larger  receiver  optics. 

A  similar  consideration  is  involved  in  evaluating  the  perfor¬ 
mance  of  an  optical  heterodyne  detector,  which  involves  a  double 
integral  of  the  MTF  over  the  receiving  aperture.  Based  on  an  MTF 
which  rapidly  tends  to  zero,  Fried^^^^  has  predicted  a  limit  to  the 
achievable  average  signal-to-noise  ratio,  no  matter  how  large  the 

The  quantity  is  the  propagation  distance  where  the  average 
field  is  down  by  the  factor  e  ^  and  is  determined  from  the  equation 
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detection  collection  aperture.  A  high  saturation  value  o£  the  MTF 
suggests  that  further  improvement  of  the  performance  with  larger 
apertures  might  still  be  a  reasonable  consideration  for  sufficiently 
short  paths. 
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